Magnetic structure of EuFe2 As2 determined by single crystal neutron diffraction 
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Among various parent compounds of iron pnictide superconductors, EuFe2As2 stands out due to the presence 
of both spin density wave of Fe and antiferromagnetic ordering (AFM) of the locaUzed Eu^^ moment. Single 
crystal neutron diffraction studies have been carried out to determine the magnetic structure of this compound 
and to investigate the coupling of two magnetic sublattices. Long range AFM ordering of Fe and Eu spins 
was observed below 190 K and 19 K, respectively. The ordering of Fe^^ moments is associated with the wave 
vector k = (1,0,1) and it takes place at the same temperature as the tetragonal to orthorhombic structural phase 
transition, which indicates the strong coupling between structural and magnetic components. The ordering of Eu 
moment is associated with the wave vector k = (0,0,1). While both Fe and Eu spins are aligned along the long 
a axis as experimentally determined, our studies suggest a weak coupling between the Fe and Eu magnetism. 

PACS numbers: 75.25.+Z, 75.50.Ee; 74.70.-b 



I. INTRODUCTION 

The recent discovery of pnictide superconductors has 
drawn extensive attention because it provides a new op- 
portunity to investigate the mechanism of superconductivity 
ill I2I ISl y, [st] . Most of the research on pnictide supercon- 
ductors has focused on /?FeAs(0;cFi-x)(with R = La, Nd and 
Sm) and AFe2As2 (with A = Ba, Ca and Sr), the so called 
' 11 1 1 ' and '112' families. These two families are closely re- 
lated since both of them adopt a layered structure with a sin- 
gle FeAs layer in the unit cell of ' 1 1 1 1 ' and two such layers 
in the unit cell of '122'. The superconducting state can be 
achieved either by electron or hole doping of the parent com- 
pounds iS H Si- Till now, the highest Tc attained is 57.4 
K in the electron doped Cao.4Nao.6FeAsF '1111' compound 
JTo'], while for '122' family the highest of 38 K is reached 
in the hole doped Bao.6Ko.4Fe2As2 [9]. Considering that the 
electronic states near the Fermi surface are dominated by con- 
tributions from Fe and As, it is believed that the FeAs layers 
are responsible for superconductivity in these compounds. 

Recent neutron diffraction experiments reveal that the for- 
mation of the spin density wave (SDW), originating from the 
long range antiferromagnetic (AFM) order of the Fe moments 
at low temperature, seems to be a common feature for all the 
iron pnictide parent compounds The on- 

set of the AFM order is also accompanied by the Tetragonal- 
Orthorhombic (T-0) structural phase transition in the '122' 
family and preceded by the T-0 phase transition for the ' 1 1 11 ' 
family. Phase diagrams of some iron pnictides clearly show 
that the magnetic order is suppressed with appropriate charge 
carrier doping of parent compound. Concomitantly, super- 



conductivity emerges and reaches a high Tc at optimal dop- 
ing ifldl . thus exhibiting features similar to high Tc cuprates 
[17]. It is generally believed that the superconductivity in iron 
pnictides is unlikely due to simple electron-phonon coupling, 
as demonstrated from extensive studies of phonon dynam- 
ics Magnetism seems to play a crucial role in the 
appearance of superconductivity and AFM spin fluctuations 
have thus been suggested to be a possible paring mechanism. 
Strong evidence on the presence of resonant spin excitation 
in the superconducting phase has indeed been obtained from 
recent inelastic neutron scattering experiments on several op- 
timally doped '122' superconductors 

EuFe2As2 is a peculiar member of iron arsenide AFe2As2 
family since the A site is occupied by Eu^^, which is an 5"- 
state (orbital angular momentum L = 0) rare-earch ion pos- 
sessing a 4/^ structure with the electron spin S = 7/2. The 
theoretical efl'ective magnetic moment of Eu^^ ion is 7.94 yUg. 
As revealed by Mossbauer and magnetic susceptibility mea- 
surement on single crystals, the Eu^^ moments order antifer- 
romagneticlly below ~ 20 K 13, [Hi- It is also reported 
that the moment of Eu^^ can be realigned ferromagneticly by 
applying a magnetic field Besides, superconductivity 

can also be achieved by replacing Eu by alkali metals, e.g. the 
Tc is observed to be 31 K and 34.7 K for Euo.5Ko.5Fe2As2 
and Euo.7Nao.3Fe2As2 [29], respectively. Unlike BaFe2As2, in 
which the superconductivity emerges with the Ni substitution 
of Fe, the SDW is suppressed in EuFe2-xNi;cAs2 without the 
emergence of superconductivity 1301 . Furthermore, the mag- 
netic ordering of Eu^^ moments evolves from AFM to ferro- 
magnetic at higher levels of Ni doping. 

Since magnetism and superconductivity appears to be inti- 
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mately related in iron pnictides, it is therefore equally impor- 
tant to understand the magnetic properties especially for the 
compounds that contain the magnetic lanthanide ions. The in- 
vestigation of the interplay between the lanthanide and iron 
magnetism may also be crucial for a deeper understanding of 
the magnetic and electronic properties of iron pnictides. For 
EuFe2As2, the magnetic ordering and the details of magnetic 
structure have not been clarified so far via single-crystal neu- 
tron diff'raction due to the extremely large neutron absorption 
cross-section of Eu. Here we report neutron diffraction stud- 
ies on a high-quality EuFe2As2 single crystal using the hot 
neutron source. It has been observed that both the Fe^^ and 
Eu^^ moments are ordered antiferromagnetically below 190 
K and 19 K, respectively. A unique propagation vector k = 
(1,0,1) is determined for the Fe magnetic sublattice with the 
moment aligned along the a axis. Furthermore, the magnetic 
propagation vector is determined to be k = (0,0, 1) for the Eu^^ 
moment, which is also aligned along the a axis. The coupling 
between the Fe and Eu magnetic sublattices has been found to 
be rather weak. The determination of the magnetic structure 
of EuFe2 As2 would pave the way for further investigations of 
EuFe2As2 under high pressure and strong magnetic fields. 



II. EXPERIMENT 

EuFe2As2 single crystals were grown by the ffux method. 
A small amount of powdered single crystal was examined by 
means of x-ray powder diff'raction (XRD) analysis. The XRD 
pattern reveals a single phase of EuFe2As2 in the tetragonal 
ThCr2Si2 structure with space group I4/mmm at room tem- 
perature. The samples have also been characterized via the 
measurements of heat capacity, resistivity and magnetic sus- 
ceptibility. Two prominent phase transitions can be identified 
respectively at 190 and 19 K, consistent to those previously 
A 50 mg single crystal with dimension 



reported 

about 5x5x1 mm^ was selected for neutron diffraction 
experiment, which was performed on hot-neutron four-circle 
diff'ractometer HEIDI at FRM II, Garching (Germany). A Cu 
(220) monochromator was selected to produce a monochro- 
matic neutron beam with the wavelength at 0.868 X. An 
Er filter was used to minimize the A/2 contamination. Sin- 
gle crystal sample was mounted on a thin aluminium holder 
inside a standard closed-cycle cryostat. The diff'raction data 
were collected using a ^He single detector at diff'erent tem- 
peratures from 300 K down to 2.5 K. A fine collimation (~ 
150 in front of the sample and a narrow opening of the de- 
tector slits were adopted to achieve a suflftcient resolution, in 
order to determine precisely the structural splitting due to or- 
thorhombic twinning and magnetic modulation wave vectors. 
To ensure the inclusion of the contributions from all possi- 
ble twinned domains, the integrated intensities were collected 
with the setup adopting a 60' collimation and an angular open- 
ing of both horizontal and vertical detector slits set at 4.5 
degree. Furthermore, the integrated intensities for the reffec- 
tions with 20 > 60° and 26 < 60° were obtained respectively 



^ 7000 

Oj 6000 
W 

§ 4000 
O 

^ 3000 

^ 2000 
c 

03 1000 



"J 


\ (220)^ 




17.7 18.0 18.3 18.6 18.9 19.2 

CO (degree) 


(C) A 

D1+D3 i \i 


I \ D2+D4 



^ 12000 

g3 10000 

"2 8000 
c 

O 6000 

^ 4000 

§ 2000 



(b) 1 


\ (200), 





12.0 12.5 13.0 

CO (degree) 



17.7 18.0 18.3 18.6 

CO (degree) 




24.5 25.0 25.5 26.0 26.5 27.0 

CO (degree) 



(e) 



a ■ 




b 














Re 






1 


























,-'(110) 








(110j\ 










D1 
















D4 






























D2 




I fab 






D3 








b 


a 











k(r. I. u.) 



FIG. 1: (Color online) (a)(b) Omega scans of tetragonal (220)r and 
(400)r nuclear reflection at 300 K, respectively. (C) Omega scan 
of orthorhombic (400)^ nuclear reflection at 2.5 K. The splitting of 
the reflection indicates the existence of twining, (d) Triple splitting 
of the rocking curve of orthorhombic {220)o reflection at 2.5 K, as 
measured with the (hkO) aligned nearly in the horizontal scattering 
plane (e) Schematic presentation of the twinned orthorhombic lattice 
in real space. Four domain patterns are marked as D1-D4. a and b de- 
note long and short axis of orthorhombic lattice. Red arrows indicate 
the direction of the Fe magnetic moment, (f) Schematic presentation 
of the reciprocal space corresponding to the twinned orthorhombic 
domains, (g) The contour map of orthorhombic (400) ^ reflection at 
2.5 K. (h) Q scan of (400)^ reflection. 



via the 6-26 and the rocking-curve scans. The obtained data 
used for the structural refinements were normalized by mon- 
itor counts and corrected for the Lorentz factor. DATAP pro- 
gram is used to carry out absorption correction by considering 
the size and shape of crystal |31]. The absorption coeflftcient 
yu is calculated to be 2.61 mm"^ and the transmission factors 
are deduced to be only in the range from 2.1% up to 14.2% 
due to the extremely strong absorption. Determination of both 
the nuclear and magnetic structures was performed by using 
the FULLPROF program suit |32]. The scale factor derived 
from the crystal structure refinement was used to determine 
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FIG. 2: (Color online) (a) The schematic diagram of (h I) plane 
in the first quadrant of reciprocal space. The circular, rhombic and 
square symbols represent the nuclear reflection as well as the mag- 
netic reflection attributed to Fe and Eu magnetic sublattices. (b) Long 
/ scan on (01/), (02/) and (11/) reflections. 



the magnitude of magnetic moment from the magnetic reflec- 
tions. 



III. RESULTS AND DISCUSSION 

First of all, the crystal structure of EuFe2As2 is described 
within the framework of tetragonal symmetry at 300 K. The 
oj scans of selected nuclear (220)r and (200)r reflections with 
mosaic width of ~0.22° are shown in Fig. 1 (a) and (b), which 
indicate the good quality and homogeneity of the single crys- 
tal. Upon cooling down, a splitting is observed for orthorhom- 
bic (400)o and (040)o reflections (Fig. 1(c)). Note that those 
two reflections are corresponding to the (220) r in tetragonal 
setting according to the Tetragonal- Orthorhombic symmetry 
relation. The observed splitting of (400)o is the indication 
of T-0 structural transition and accompanied twinning con- 
figuration due to the interchange of the orthorhombic a and 
b axes. It is known that twinning in orthorhombic structure 
will result in four diflTerent domain patterns (SSlls?], as illus- 
trated in Fig. 1 (e). Two of domains shared the same (110) 



plane and formed the domain pairs (Dl and D2), while an- 
other two shared the (110) plane (D3 and D4). In principle, it 
is possible to observe single peak, two or three or four peaks 
depending on the selected reflections and the resolution of the 
instrument. In Fig. 1(c), the left and right peaks can be as- 
signed to the contributions from the domains (D1-I-D3) and 
(D2-hD4) respectively. Note that the 6i;-scan is performed with 
open detector slits. Two Gaussian peaks were used to fit the 
(400) (9 and (040) o reflections and the domain population ra- 
tio is estimated to be around 1 : 1 for the (hOO) and (OkO) twins, 
i.e. Dl-hD3 ^ D2-hD4. The a;-scan of (220) is examined after- 
ward with the (hkO) aligned in the horizontal scattering plane 
to obtain more detailed information about domain population 
(Fig. 1(d)). The occurrence of twinning and T-0 structural 
phase transition can be confirmed from the clear presence of 
triple splitting of (220) nuclear reflection. Usually, the reflec- 
tions with h = k ^ are triply split in twinned orthorhombic 
lattice and the peak in the center is attributed to pairs that share 
the same (T 10) or (1 10) plane, while the peak at left and right 
sides corresponds to the rest two domains. In Fig. 1(d), it 
can be seen that all three peaks showed almost equal intensity. 
This strongly indicates that the domain population exists the 
following relationship: Dl ^ D3+D4 ^ D2. Hence the do- 
main population for all those four diflTerent domain patterns 
can be determined roughly as 2:2:1:1. In order to investigate 
the distribution of nuclear reflection in reciprocal space and 
determine the lattice parameter accurately, two dimensional 
plot of (400) o reflection is shown in Fig. 1 (g). The splitting 
of (400)(9 can also be clearly seen. Totally 280 nuclear reflec- 
tions were collected for nuclear structure refinement within 
the Fmmm space group. Several strong reflections were ex- 
cluded from the refinement because of the significant extinc- 
tion. All atoms were refined with the isotropic temperature 
factor. The refinement results of crystal structure are listed in 
Table 1. The lattice parameters are deduced tohea = 5.537(2) 
\,b = 5.505(2) X and c = 12.057(2) X at2.5 K, which are 
in good agreement with a previous report [35]. 

To clarify the magnetic structure of EuFe2 As2 at low tem- 
perature, the sample was cooled to 2.5 K, which is well be- 
low the reported Fe^^ and Eu^^ magnetic ordering temper- 
atures. Considering the existence of the twined (hOO) and 
(0^0) domains, extensive search of magnetic reflections was 
performed in the a^-c^ reciprocal space as schematically il- 
lustrated in Fig. 2(a). Additional search was also performed 
in the (hhl) reciprocal plane. Fig. 2(b) shows three typical 
long / scans in the reciprocal space where in addition to the 
expected nuclear reflections, two sets of magnetic superstruc- 
ture reflections can be clearly identified with two magnetic 
propagation wave vectors (1,0,1) and (0,0,1) respectively. As 
an example, Q scan of (101)m magnetic reflection is plotted 
in Fig. 3(a) and the same scan at 200 K is also plotted to- 
gether for comparison. In Fig. 3(b), the 0-26 scan of nuclear 
(002)a^ and magnetic (003 )m reflections show the same peak 
center, which indicates that the magnetic structure is commen- 
surate in nature. The contour map of (103)m and (401)m re- 
flections fully illustrated the intensity distribution as shown 
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FIG. 3: (Color online) (a) The comparison of the Q scan of (101)m 
magnetic reflection at 2.5 and 200 K. The (101)m reflection is ob- 
served in k scan because of the existence of twining, (b) The 6-26 
scan of (003)a^ nuclear and (003)m magnetic reflections at 2.5 K, the 
same scan of (003)m magnetic reflection at 300 K is also plotted for 
comparison, (c) The contour map shows the Q dependence of the 
(103)m magnetic reflection, (d) Q scan of (103)m magnetic reflec- 
tion. The (103)m reflection is observed in k scan because of the exis- 
tence of twining, (e) The contour map of (041)^ magnetic reflection 
indicates the contribution of the magnetic reflection of Eu magnetic 
sublattice. (f) Q scan of (04 1)^ magnetic reflection. 

in Fig. 3(c) and Fig. 3(e). As already discussed, the contour 
map of (400)(9 nuclear reflection (Fig. 1(g)) clearly shows two 
reflections attributed to the (/zOO) and (0^0) twins. Two peak 
centers with k = 3.967 and 4.01 can be obtained by fitting the 
Q scan of (400)o reflection(Fig. 1(h)). In Fig. 3.(d), the Q 
scan of (103)m reflection can be fitted by a single Gaussian 
function with k = 0.991. This strongly indicates that the (hOl) 
type reflections (with h and / equal to odd number) are asso- 
ciated with the (hOO) domain and they can thus be described 
with the propagation wave vector k = (1,0,1). This wave vec- 
tor is exactly the same as observed in other '122' pnictides, 
such as BaFe2As2 l[l3ll and CaFe2As2 which is related 
to the AFM order of Fe^^ moments. The magnetic structure 
refinement was then carried out to determine the magnitude 
and direction of Fe^^ moment. The magnetic structure with 
Fe saturation moment of 0.98(8) yUg aligned along the long a 
axis is deduced. Note that the origin of AFM order in FeAs- 
based pnictides is still a matter of controversy. It is argued 
that the AFM order of Fe^^ may arise from the SDW order 
due to the Fermi surface instability under the itinerant model 



FIG. 4: (Color online) (a)-(e) Omega scans of series of (02/) (with / = 
odd) reflections at 2.5 K. The integrated intensities of (20/) and (02/) 
can be obtained by fitting the curves with two Gaussian functions, (f) 
The ratio between (20/) and (02/) reflections shows good agreement 
with the calculated value. 



TABLE I: Refined results of the crystal and magnetic structures for 
EuFe2As2 at 2.5 K (space group: Fmmm, Z = 4). 



Atom/site 


X 


y 


z 




Eu (4a) 











0.81(3) 


KMaiMB) 


(0,0,1), 6.8(3) 








Fe (8/) 


0.25 


0.25 


0.25 


0.26(3) 


k, MaiMB) 


(1,0,1), 0.98(8) 








As (8/) 








0.363(5) 


0.25(3) 



a, b,c(\): 5.537(2), 5.505(2), 12.057(2) 

Number of reflections (Nuclear): 280 

RF^, RF^^, RF(%),x^ Nuclear: 9.34, 9.67, 6.22, 7.1 

Number of reflections (Magnetic): 228 

RF^, RF^^, RF(%lx^ Magnetic: 9.42, 7.68, 6.53, 5.7 



fsS] . While other evidences support that the AFM order has 
a local moment origin as in Mott insulator [38|], such as 
the parent compound of high Tc cuprates. Recent spin wave 
excitation study on CaFe2As2 suggests that the magnetism of 
iron arsenide might be resulted from a complicated mixture of 
localized and itinerant properties and it should be understood 
by considering both the localized and itinerant electrons | 39]. 

Consequently, the magnetic reflections with a propagation 
wave vector k = (0,0,1) (with h even and / odd) are due to the 
long range order of the localized Eu^^ moments. However, 
the moment direction of Eu^^ moments is still indeterminate. 
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FIG. 5: (Color online) Integrated intensities of the nuclear and mag- 
netic Bragg reflections collected at 2.5 K are plotted against the cal- 
culated values. See text for details of the crystal and magnetic struc- 
ture models. 



Symmetry analysis based on the representation theory indi- 
cates that the magnetic representation T for magnetic Eu^^ on 
4a site is decomposed into three one dimensional irreducible 
representations: Fi, F2 and F3. The Eu^^ moments are aligned 
in the c, a or b direction according to those three representa- 
tions. The observation of nonzero intensity of (00/)(with / = 
odd) reflections clearly exclude the representation Fi . The oj 
scans on several (hkO)M (with both h and k = odd) reflections 
also exhibit considerable intensity. Thus, the moment of Eu^^ 
is expected to be aligned either along the a orb direction in the 
ab plane. The Q scan on (041)m reflection (Fig. 3(f)) giving a 
peak position of ^ = 4.01, which is exactly equal to the larger 
k value of (400)0 nuclear reflection. Therefore, the moment 
direction of Eu^^ can be determined as along the a direction 
since the intensity ratio between (041)m and (401)m magnetic 
reflections approximate equals to 73:1 for this arrangement. 
The structure mode is confirmed further by oj scan of series of 
(02/) (with / = odd) reflections as shown in Fig. 4. Similar to 
some nuclear reflections, both (20/) and (02/) magnetic reflec- 
tions was detected due to the twinning configuration. How- 
ever, the intensity ratio between (20/) and (02/) changes grad- 
ually with the change of the angle between the scattering plane 
and the c axis. The calculated intensity ratio of (20/)/(02/) 
for diff'erent / are plotted in Fig. 4(f) and it agrees well with 
the observed values which derived from the oj scans directly. 
By taking into account of twinning components properly, the 
refinement on Eu^^ magnetic sublattice was carried out with 
the aforementioned magnetic structure model. The calculated 
structure factors are plotted against those observed and shown 
in Fig. 5. The reliable agreement factors confirms the pro- 
posed magnetic structural model eventually, i.e. the Eu^^ mo- 
ment aligns along a direction with the wave vector k = (0,0,1) 
and magnitude of 6.8(3) yUg. Thus the magnetic structure of 
EuFe2As2 is unambiguously determined as illustrated in Fig. 
6. The moment direction of Eu^^ is also consistent with our 
resonant x-ray scattering (RXS) measurement [40]. 

Fig. 7(a) shows the temperature dependence of the (112)m 




^ As 



FIG. 6: (Color online) Illustration of the magnetic structures of 
EuFe2As2 at 2.5 K. The Fe moments align along a direction and order 
antiferromagnetically in both a and c directions. The Eu moments 
align along a direction and order antiferromagnetically in c direction 
only. The gray line outlines the orthorhombic unit cell. 



and (003 )m magnetic reflections attributed to the ordering of 
Eu^^ moments. The onset temperature of Eu^^ magnetic order 
is deduced to be 19 K, which is in good agreement with previ- 
ous report on electronic and magnetic measurements | 24, ■2?]. 
The magnetic ordering temperature of Fe^^ moment is esti- 
mated to be 190 K based on the temperature dependence of 
the (101)m and (103)m magnetic reflections (see Fig. 7(b)). 
The T-0 structural phase transition also takes place at 190 K 
as revealed by the sharp change of full width at half maximum 
(FWHM) in (040) o nuclear reflection. First principle calcula- 
tions suggest that the nearest and next nearest neighbor su- 
perexchange interactions between Fe ions lead to a frustrated 
magnetic ground state in pnictides with parallel and antipar- 
allel arrangement of Fe spins in FeAs layer ll4l|] . Usually, 
the magnetic frustration can be lifted by a structural distortion 
from low symmetry to high symmetry phase. This may be the 
origin of the strong coupling between the structural and mag- 
netic phase transitions observed in EuFe2As2 and other iron 
pnictides lfT2[ [Bl H . 

Due to the localized nature of Eu 4/ state, the AFM cou- 
pling of Eu^^ moments would be described by the indirect ex- 
change, e.g. the Ruderman-Kittel-Kasuya-Yosida(RKKY) in- 
teraction as suggested by Ren et al [30]. Besides of the Eu-Eu 
and Fe-Fe interactions, the strength of the interaction between 
the Eu and Fe magnetic sublattices is also an interesting issue. 
Similar to CaFe2As2, the SDW transition in EuFe2As2 can 
also be suppressed continuously by applying the pressure due 
to the weakening of nearest Fe-Fe exchange coupling [42]. 
Whereas the AFM ordering temperature of Eu sublattice does 
not change significantly even the compound exhibits the pos- 
sible reentrant superconducting state. This may suggest the 
weak interaction between the Eu and Fe magnetic sublattices. 
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FIG. 7: (Color online) (a) Temperature dependence of integrated in- 
tensity of (040)o nuclear reflection as well as the (103)m , (112)m and 
(003) magnetic reflections below 22 K. (b) Temperature dependence 
of integrated intensity of (103)m and (101)m reflections; temperature 
dependence of FWHM of (040)^ reflection. 

which is supported by the full potential electronic structure 
calculation [43] . In present neutron work, we also did not ob- 
serve any detectable change in Fe^^ magnetic moment when 
temperature passing through the Eu^^ magnetic ordering tem- 
perature (Fig. 7(a)). The coupling between Fe^^ and Eu^^ 
moments in the ordered state varnishes entirely within a mean 
field model due to geometrical frustration. Parallel alignment 
of the Fe^^ and Eu^^ moments might be resulted from thermal 
or ground state fluctuations, as suggested in an order by dis- 
order scheme 13, Es]]. Those results are in contrast to some 
' 11 11 ' compounds, such as PrFeAsO, in which the interplay 
between Fe and Pr ordering moments might drive the negative 
thermal expansion [46] . 

IV. CONCLUSION 

Single crystal neutron difl'raction experiment using a hot 
neutron source was performed to investigate the crystal and 
magnetic structure of EuFe2As2. With decreasing temper- 
ature, the antiferromagnetic order of Fe^^ moments set in 
at 190 K with the propagation vector k = (1,0,1). Similar 
to BaFe2As2 and CaFe2As2, the tetragonal to orthorhombic 
structural transition occurs simultaneously with the AFM or- 
der, which indicates the strong coupling between the lattice 
and Fe magnetic degree of freedom. Below 19 K, the Eu^^ 
moments order antiferromagneticly with the propagation vec- 
tor k = (1,0,1) and are aligned along the a axis, our studies 
also suggest a weak coupling between the Fe^^ and Eu^^ mag- 
netic sublattice. 
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